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PREFACE 

The  study  reported  herein  was  performed  under  Department  of  the  Army 
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the  Army  (R&D). 
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Chief,  Rock  Mechanics  Applications  Group  (RMAG),  Engineering  Geology  and  Rock 
Mechanics  Division  (EGRMD),  Geotechnical  Laboratory  (GL),  and  un^ler  the 
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III  was  Chief,  GL.  Ms.  M.  E.  Glynn,  RMAG,  assisted  in  the  collection  and 
reduction  of  test  data.  Ms.  M.  A.  Kirklin,  EGRMD,  typed  the  initial  draft 
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COL  Til  ford  C.  Creel,  CE,  and  COL  Robert  C.  Lee,  CE,  were  Commanders 
and  Directors  of  WES  during  the  conduct  of  the  study.  COL  Allen  F.  Grum, 

USA,  was  Director  of  WES  during  the  preparation  and  publication  of  this  re¬ 
port.  Mr.  Fred  R.  Brown  and  Dr.  Robert  W.  Whalin  were  Technical  Directors. 


% 


r,- 

C** 

tV 


V  --J.  -  J 


CONTENTS 


Page 


PREFACE  .  1 

PART  I:  INTRODUCTION .  3 

Background  .  3 

Scope .  5 

PART  II:  CONSTITUTIVE  RELATIONSHIPS  .  6 

Current  Practice  .  6 

Proposed  Constitutive  Model  .  10 

PART  III:  PROPOSED  STRAIN  BASED  FAILURE  CRITERION  .  15 

Background .  15 

Triaxial  Test  Data .  18 

Correlations  Between  Maximum  Shear  Stress  and  Axial  Strain.  .  19 

Correlations  Between  Minor  Principal  Stress  (Confining 

Stress)  and  Axial  Strain  .  24 

Deformation  Modulus  .  27 

PART  IV:  SUMMARY  AND  CONCLUSIONS .  30 

Summary .  30 

Conclusions .  30 

REFERENCES .  32 

TABLES  1-3 

APPENDIX  A:  SANDSTONE  AND  GRANITE  TEST  DATA .  A1 


2 


ROCK  FAILURE  CRITERIA  FOR  DESIGN  OF 
UNDERGROUND  STRUCTURAL  SUPPORTS 


PART  I:  INTRODUCTION 
Background 

Failure  Criteria 

Bieniawski  (1984)  defines  a  criterion  of  failure  as  "an  algebraic  expres¬ 
sion  of  the  mechanical  condition  under  which  a  material  fails  by  fracturing  or 
deforming  beyond  some  specified  limit.  This  specification  can  be  in  terms  of 
load,  deformation,  stress,  strain,  or  other  parameters."  Bieniawski  listed  no 
less  than  fourteen  failure  criteria  that  have  been  used  in  the  past  or  are 
currently  being  used  to  define  failure  for  rock.  These  criteria  are  listed  in 
Table  1.  While  the  empirical  rock  mass  strength  criterion  proposed  by  Hoek  and 
Brown  (1980)  (see  Table  1)  attempts  to  define  failure  in  rock  masses,  the  major¬ 
ity  of  the  criteria  listed  in  Table  1  were  either  specifically  developed  for 
i  itoct  rock  or  adopted  to  define  intact  rock  failure. 

Table  1  was  not  intended  to  be  all  inclusive.  In  addition  to  the  cri¬ 
teria  listed,  criteria  are  available  for  a  third  potential  mode  of  failure  - 
sliding  along  a  single  discontinuity.  The  better  known  of  the  single  plane 
sliding  failure  criteria  include  the  bilinear  criterion  proposed  by  Patton 
(1966)  and  Goldstein  et  al.  (1966)  and  the  curvilinear  criterion  proposed  by 
Barton  (1971  and  1973),  and  Ladanyi  and  Archambault  (1969). 

The  mechanical  complexity  of  failure  increases  from  intact  rock  to  dis¬ 
continuous  rock  to  rock  mass.  Although  much  progress  has  been  made  over  the 
past  200  plus  years  since  Coulomb  proposed  his  criterion  for  failure  in  1773, 
the  practicing  engineer  is  still  in  need  of  a  general  criterion  that  better 
defines  rock  failure.  However,  because  of  basic  differences  in  mechanical 
behavior  between  the  potential  modes  of  rock  failure,  it  is  doubtful  that  a 
single  theoretical  criterion  will  ever  be  developed  to  describe  all  three  modes 
of  failure.  Bieniawski  (1984)  is  of  the  opinion  that  to  meet  the  immediate 
needs  of  the  practical  rock  engineer,  attention  should  be  directed  to  empiri¬ 
cal  criteria  for  estimating  the  triaxial  strength  of  rock.  Such  criteria  can 
be  selected  by  fitting  a  suitable  equation  into  experimental  dat?  and  they 
need  not  have  a  theoretical  basis. 

Because  of  the  complexities  associated  with  rock  failure,  it  is  doubtful 


that  failure  criteria  founded  on  theoretical  considerations  will  be  forth¬ 
coming  in  the  near  future.  In  this  respect,  empirical  criteria  offer  a 
viable  interim  approach  for  predicting  stress-strain  responses  of  deep  under¬ 
ground  excavations  such  as  the  Deep  Basing  concept  being  considered  for 
deployment  of  the  MX  missile. 

Existing  Strain  Based  Criteria 

An  extensive  literature  review  located  only  one  failure  criterion  based 
on  strain.  The  criterion  was  proposed  by  St.  Verrant  in  1870  (Nadai,  1950) 
and  became  known  as  the  Maximum  Elastic  Strain  Theory  or  the  Equivalent  Stress 
Theory.  The  criterion  was  developed  for  expressing  the  condition  of  brittle 
failure  and  even  yield  for  metals.  According  to  this  theory,  the  maximum 
positive  elastic  extension  of  the  material  determines  failure  of  either  kind. 

From  Hooke's  Laws  of  Elasticity,  if  the  elastic  strain  defined  by  Equation  1 
is  positive,  (sign  convention  tension  is  positive) 

=  r  [oi  -  (o2  +  03)3  £  (1) 

u 

where  ei  =  Major  principal  strain 
E  =  Modulus  of  elasticity 

oi,  02,  03  =  Major,  intermediate,  and  minor  principal  stresses,  respectively. 
V  =  Poisson's  ratio 

then  the  expression 

oi  -  V  (02  +  03)  <  Oq  (2) 

Should  be  smaller  than  the  positive  limiting  "equivalent"  (fracture  or  yield) 
stress  CTq  . 

Although  the  Maximum  Elastic  Strain  Theory  enjoyed  considerable  attention, 
particularly  in  France  and  Germany,  around  the  turn  of  the  century,  the 
theory  as  a  failure  criteria  is  not  suited  for  rock.  Consider  the  case  of 
simple  uniaxial  compression  (i.e.  02=03=  0,  -03)  where  from  Hooke's  Laws 

the  lateral  strain  is  positive  according  to: 

1  V  Oi 

£2  =  ^  [  -  v(-Ol)]  =  -  (3) 

where  £2  =  Intermediate  principal  strain. 

According  to  this  criterion  for  uniaxial  compression  v  03  <  o^  or 
oi  <  %  (from  Equations  1  &  2  for  pure  tension  03  ^  oq).  Hence  it  can  be  seen 

V 
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that  the  equivalent  fracture  stress  in  compression  is  greater  by  a  factor  of 
oq/v  than  the  fracture  stress  in  tension.  Laboratory  experience  with  rock 
indicates  that  a  typical  value  for  Poisson's  ratio  is  on  the  order  of  0.20.  In 
such  a  case,  the  criterion  would  predict  that  the  uniaxial  compressive  strength 
is  5  times  greater  than  the  tensile  strength.  However,  laboratory  test  data 
indicate  chat  the  uniaxial  compressive  strength  for  most  intact  rock  typi¬ 
cally  varies  from  10  to  15  times  the  tensile  strength. 

The  primary  short  comings  of  the  Maximum  Elastic  Strain  Theory,  as  applied 
to  rock,  is  that  the  Theory  assumes  that  the  material  behaves  in  a  perfectly 
elastic  manner.  Rock,  in  particular  a  rock  mass,  can  seldom  be  expected  to 
behave  as  an  elastic  material. 

Scope 

This  report  briefly  examines  the  role  of  constitutive  stress-strain 
relationships  as  necessary  input  to  numerical  analysis  of  rock  wall/support 
responses  of  underground  excavations.  Currently  used  models,  as  well  as  a 
proposed  constitutive  model,  are  discussed  in  terms  of  their  applicability 
and  limitations.  The  proposed  constitutive  model  requires  the  development 
of  a  failure  criterion  that  defines  strain  at  failure.  Toward  this  end,  the 
feasibility  of  developing  a  strain  based  failure  criterion  for  two  intact  rock 
types  are  investigated.  The  two  rock  types  are  intact  sandstone  and  intact 
granite.  Data  for  the  two  rock  types  were  obtained  from  triaxial  tests  results 
published  in  the  literature. 
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PART  II:  CONSTITUTIVE  RELATIONSHIPS 


Current  Practice 


For  convenience  of  discussion,  one  can  divide  the  design  of  structures 
constructed  on  or  within  a  rock  mass  into  two  general  categories:  surface  or 
near  surface  structures  and  underground  excavations.  The  distinction  between 
the  two  categories  lie  in  the  manner  in  which  acting  stresses  are  generated. 
Acting  stresses  include  stresses  tending  to  cause  instability  (driving  stresses) 
and  stresses  acting  normal  to  the  potential  failure  surface  (normal  stresses). 
For  surface  or  near  surface  structures  (e.g.  gravity  dams  and  slopes)  acting 
stresses  are,  as  a  rule,  independent  of  strain.  Resistance  to  the  stresses 
acting  on  the  structure  is  generated  by  the  inherit  strength  of  the  rock  mass 
or  interface  between  the  rock  mass  and  the  structure.  While  resisting  stresses 
in  such  cases  are  dependent  upon  strain,  satisfactory  stability  determinations 
can  be  obtained  by  ignoring  this  dependency  upon  strain  and  specifying  in  the 
stability  analyses  the  appropriate  upper  limit  of  resisting  stress  (shear 
strength)  that  can  be  developed.  Upper  limits  of  resisting  shear  stress  are 
usually  defined  by  a  failure  criterion  which  relates  either  the  principle 
stresses  or  the  normal  stresses  to  the  shear  strength  that  can  be  developed. 
Suitable  factors  of  safety,  based  on  experience,  are  specified  for  the 
solution  of  stability.  Factors  of  safety  primarily  account  for  the  uncer¬ 
tainty  of  the  assigned  upper  limit  of  resisting  shear  strength.  However, 
for  structures  that  can  tolerate  only  small  deformations  during  operation, 
appropriate  factors  of  safety  also  minimize  the  amount  of  strain  necessary  to 
develop  the  resisting  stress  required  by  the  design. 

In  underground  excavations,  as  the  preexisting  rock  support  is  removed 
during  excavation,  displacements  occur  in  the  rock  surrounding  the  opening, 
and  stresses  are  redistributed.  Hence  stress  acting  within  the  wall  rock  and 
on  the  structural  supports  are,  in  part,  a  function  of  the  relaxation  (strain) 
process.  Mechanisms  of  interaction  between  the  rock  walls  and  supports  are 
complex.  Because  of  the  complexities  the  design  processes  for  structural 
support  rely  almost  entirely  on  empirical  rock  mass  classification  systems 
(i.e.  Terzaghi  (1946),  Bieniawski  (1973),  Barton,  Lien,  and  Lunde  (1974)  and 
others).  This  does  not  mean  to  imply  that  numerical  approaches  are  not 
available  for  modeling  rock  wall/support  responses. 

Numerical  approaches  can  be  separated  into  two  general  groups  ■  continuum 
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and  discontinuum.  Continuum  approaches  include  the  finite  element,  finite 
difference,  and  boundary  element  methods.  Discontinuum  models  feature 
numerical  approaches  involving  equation  of  motion  of  particles  or  blocks. 

Such  models  are  commonly  referred  to  as  discrete  element  methods.  A  con¬ 
tinuum  approach  treats  the  rock  mass  as  a  continuum  intersected  (in  some 
models)  by  a  number  of  discontinuities,  while  discontinuum  approaches  treat 
the  rock  mass  as  an  assemblage  of  independent  blocks.  Discontinuum  approaches 
are  primarily  used  when  analyzing  the  displacement  of  independent  and  recog¬ 
nizable  rock  blocks.  Of  the  numerical  approaches  available,  the  finite 
element  method  is,  by  far,  the  most  popular  and  perhaps  the  most  potentially 
useful  analytical  design  tool. 

The  finite  element  concept  as  applied  to  continuum  mechanics  was  intro¬ 
duced  by  Turner  et  al .  (1956)  for  the  analysis  of  aircraft  frames.  The  first 
practical  application  of  the  method  to  geotechnical  problems  known  to  the 
author  was  made  by  Clough  (1962)  for  the  Little  Rock  District,  U.  S.  Army 
Corps  of  Engineers,  in  a  study  of  Norfolk  Dam.  Since  the  early  days  of  develop¬ 
ment,  the  finite  element  method  has  evolved  to  become  quite  sophisticated. 

Numerical  modeling  of  the  relaxation  process  dictates  the  use  of  con¬ 
stitutive  models  which  define  stress-strain  behavior.  The  U.  S.  National 
Committee  for  Rock  Mechanics  (1981)  listed  no  less  than  15  numerical  computer 
codes  suitable  for  application  to  underground  excavations.  The  codes,  con¬ 
sidered  representative  of  the  state  of  the  art,  are  capable  of  modeling  one 
or  more  of  the  following  constitutive  relationships. 

1 )  Linear  elastic 

2)  Nonlinear  elastic 

3)  Linear  viscoelastic 

4)  Nonlinear  viscoelastic 

5)  Elasto-plastic 

6)  Elasto-visco-plastic 

7)  Dilatancy 

Brady  and  St.John  (1982)  are  rightly  of  the  opinion  that  much  more  model 
application  and  verification  are  required  before  any  of  the  above  constitutive 
models  can  be  applied  for  practical  design  of  underground  excavations.  The 
reason  the  method  is  not  considered  a  reliable  design  tool  at  the  present 
does  not  lie  in  the  method's  ability  to  model  the  complexities  associated 
with  the  relaxation  process  but  rather  reflects  the  inability  of  the  geo¬ 
technical  engineer  to  provide  realistic  input  parameters. 
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The  difficulties  of  providing  realistic  input  parameters  can  best  be 
illustrated  by  considering  the  simp‘'e  linear  elasto-plastic  constitutive  model 
illustrated  in  Figure  1.  Several  alternatives  are  available  for  defining  the 
model.  A  frequently  used  alternative  specifies  the  deformation  moduli,  defin¬ 
ed  in  Equation  4,  as  input  parameters. 


where  =  Deformation  modulus 

0.  =  Deviator  stress  at  failure 

‘'f 

=  Axial  strain  at  failure  (major  principal  strain) 

Deformation  moduli  define  the  initial  slopes  of  the  stress-strain  lines.  A 
suitable  stress  based  failure  criterion  is  then  specified  to  define  the  upper 
limits  of  deviator  stress  at  failure.  A  second  alternative  specifies  both 
deviator  stress  and  strain  at  failure. 

Sood  approximations  of  input  parameters  for  either  alternative  can  be 
obtained  from  triaxial  tests  on  small  intact  or  broken  rock  specimens. 

However,  because  of  scale  effects,  stress-strain  responses  of  small  test 
specimens  seldom  model  prototype  rock  mass  behavior.  Hence,  input  parameters 
for  practical  underground  excavation  problems  involve  considerable  speculation. 
Generally,  in  practical  problems,  a  deformation  modulus  is  estimated  from  in 
situ  tests  such  as  the  plate  bearing  test,  and  the  upper  limits  of  deviator 
stress  are  defined  by  a  failure  criterion.  Such  an  approach  incorporates 
two  rather  severe  restrictions  into  the  analysis. 

First,  in  situ  tests  available  for  estimating  the  deformation  modulus 
for  rock  mass  are  typically  limited  to  an  applied  stress  range  of  0  to  6.9  MPa 
(0  to  1000  psi).  The  maximum  capacity  of  available  loading  systems  generally 
limit  the  upper  end  stress  range.  Complete  stress-strain  responses  for  most 
rock  mass  require  stress  ranges  on  the  order  of  one  magnitude  greater  than 
maximum  stresses  available  for  in  situ  testing.  In  addition,  stress-strain 
responses  are  typically  strangely  curvilinear  at  low-stress  levels  as 
illustrated  in  Figure  2.  The  strong  initial  curvature  is  primarily  due  to 
closure  of  discontinuities  (joints)  within  the  rock  mass  and  to  some  extent 
closure  of  microcracks  within  the  intact  segments  of  rock.  The  number, 
orientation,  tightness,  and  type  of  infilling  of  the  joints  control  the 
extent  of  curvature.  Modulus  values  determined  from  in  situ  tests  at 
relatively  low  stress  levels,  for  example  the  slope  of  line  OA  in  Figure  2. 


are  dilatant  at  stress  levels  approaching  failure,  at  least  for  brittle 
failure.  Hence,  one  would  expect  lower  peak  deviator  stresses  from  drained 
tests  than  from  undrained  tests.  The  fact  that  this  trend  is  not  evident 
in  Figure  6  suggests  that  the  drained  tests  specimens  might  not  have  been 
completely  saturated  or  that  reduction  in  deviator  stress  was  insufficient  to 
generate  a  noticable  trend  within  the  overall  scatter  of  data. 

Brittle-ductile  transition.  Insufficient  detailed  data  prohibit  a 
conclusive  distinction  between  brittle  and  ductile  failure.  However,  the 
limited  data  available  in  the  form  of  stress-strain  curves  suggest  that  the 
majority  of  sandstones  included  in  this  study  behaved  in  a  brittle  manner  for 
normalized  maximum  shear  stress  ratios  below  approximately  3.5.  This  maximum 
shear  stress  ratio  corresponds  to  a  normalized  axial  strain  ratio  of  approxi¬ 
mately  5.0.  A  visual  inspection  of  Figure  6  indicates  that,  over  the  apparent 
brittle  failure  range,  a  linear  correlation  would  fit  the  data  rather  well. 


Granite  Data 

Correlation.  Figure  7  shows  a  plot  of  normalized  axial  strain, 
vs  normalized  maximum  shear  stress,  (oi  -  03)720^^  ,  for  the  granite  data 
given  in  Table  A1  of  Appendix  A.  The  expression  for  the  line  best  fitting 
all  the  data  is  as  follows: 


0.26  +  0.73 


+ 


(20) 


Data  scatter.  Like  the  sandstone  data,  data  scatter  for  the  granite 
tests  was  reasonable  considering  the  sources  of  data  and  range  of  representa¬ 
tive  unconfined  compressive  strengths  (80  MPa  (11  ,600  psi)  ^  248  MPa 

(36,000  psi)).  It  is  interesting  to  note  that  the  Stone  Mountain  granite 
(and  perhaps  the  Quartz  Diorite  data)  data  in  Figure  7  deviate  from  the 
general  trend.  One  possible  explanation  for  this  deviation  might  be  that 
the  loading  frame  used  in  these  tests  may  have  been  insufficiently  stiff  for 
a  high  modulus  rock  like  granite.  A  soft  loading  frame  could  cause  dumping 
of  strain  energy  stored  in  the  loading  frame  resulting  in  a  low  observed 
deviator  stress  at  failure. 

Brittle-ductile  transition.  An  inspection  of  available  stress-strain 
curves  for  the  granite  data  given  in  Table  A1  indicated  that  all  tests  be¬ 
haved  in  a  brittle  manner.  It  is  interesting  to  note  the  almost  linear 
nature  of  Equation  20. 
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Axial  Straln/Unconfined  Compressive  Axial  Strain  -  e/e 


Since  o,  for  the  unconfined  compression  test  the  unconfined  compressive 
stress  is  indeed  the  major  principal  stress  ai  (i.e.  when  03=0  then 
oi  =  for  pure  compression  loading  (i.e.  no  tensile  stresses  generated), 
the  minimum  value  the  normalized  maximum  shear  stress  component  in  Equation  15 
can  assume  is: 


Ol  0 


0.5 


Likewise,  when  03  =  0  and  e  =  the  minimum  value  of  the  dependent 
normalized  strain  component  is: 


(16) 


^  .  1,0  (17) 

^uc  ^uc 

Equations  obtained  from  Curve  Fit  required  adjustments  so  that  minimum 
values  were  satisfied.  Putting  Equations  16  and  17  into  Equation  15  resulted 
in  required  expressions  that  assigned  values  of  constants  a,  b  and  c  must 
satisfy  for  the  unconfined  case  as  follows: 


1 .0  =  a  +  0.5b  +  0.25c  (18) 

Curve  Fit  values  of  a  ,  b  ,  and  c  were  then  adjusted  so  as  to  satisfy 
Equation  18  and  provide  a  visual  best  fit  of  the  data. 

Sandstone  Data 

Correlation.  Figure  6  shows  a  plot  of  normalized  axial  strain, 
vs  lormalized  maximum  shear  stress,  (oj  -  the  sandstone  data 

given  in  Table  Al  of  Appendix  A.  The  expression  for  the  line  in  Figure  6 
is  as  follows: 


Forcing  the  curve  to  satisfy  the  minimum  requirements  for  the  maximum  shear 
stress  and  strain  components  ((ai  -  03)720  and  e/e,,_  ,  respectively)  as 
discussed  above  shifted  the  curve  slightly  to  the  right  as  evident  in  Figure  6. 

Data  scatter  Considering  the  sources  of  data  and  the  range  of  strengths 
as  reflected  by  the  unconfined  compressive  strengths  listed  in  Figure  6  (30  MPa 
(4280  psi)  1  260  MPa  (37,500  psi)),  data  scatter  is  not  thought  to  be 

excessive.  It  should  also  be  noted  that  the  sandstone  data  include  tests  on 
wet  and  dry  specimens  as  well  as  drained  and  undrained  tests.  Robinson  (1955), 
Byerlee  and  Brace  (1967),  Dyerlee  (1975)  and  others  have  demonstrated  the 
validity  of  effective  stress  in  porous  rocks  such  as  sandstone.  Most  rocks 
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the  literature  by  projecting  from  the  point  of  failure  to  the  appropriate  axis 
defining  stress  and  strain  values.  Values  assigned  for  stress  and  strain  are 
subject  to  some  error  due  to  the  inaccuracies  of  scaling  the  data.  The  data 
represent  results  from  a  variety  of  different  testing  devices  and,  no  doubt, 
testing  procedures.  Nicholson  (1983)  is  of  the  opinion  that  test  results  can 
vary  significantly  for  tests  conducted  in  different  triaxial  devices  even 
though  care  is  taken  to  insure  that  test  procedures,  specimen  preparation,  and 
specimens  are  identical  in  so  far  as  possible.  Hence,  considering  the  variety 
of  data  sources,  one  would  expect  a  certain  amount  of  data  scatter. 

During  the  course  of  this  study  two  correlations  between  stress  and 
strain  were  found  to  be  of  potential  use  for  analyzing  underground  excavation 
behavior.  In  particular,  these  correlations  were  between  the  maximun  shear 
stress  ((oi  -  03)72)  and  axial  strain  (e)  and  between  the  minor  principal  stress 
(confining  stress  03)  and  axial  strain.  Both  stress  parameters  for  each  rock 
type  were  normalized  with  respect  to  the  unconfined  compressive  strength  for 
that  rock  type.  Likewise,  axial  strains  for  each  rock  type  were  normalized 
with  respect  to  the  axial  strain  of  the  unconfined  compressive  strength  for 
that  rock  type.  Normalization  allowed  a  relative  comparison  between  stress- 
strain  responses  of  specimen  of  the  same  rock  type  (i.e.  sandstone  and  granite) 
but  different  stress-strain  behavior.  Table  A1  (Appendix  A)  lists  the  test 
data  according  to  rock  type,  principal  stresses,  axial  strain,  and  symbols 
used  in  following  correlation  plots. 

Correlations  Between  Maximum 

Shear  Stress  and  Axial  Strain 


Curve  Fitting 

Considering  the  propensity  for  data  scatter,  a  reasonable  correlation  was 

found  to  exist  between  normalized  axial  strain  and  normalized  maximum  shear 

stress  for  both  the  sandstone  and  granite  data.  Lines  best  fitting  the  data 

were  obtained  from  a  curve  fitting  routine,  entitled  Curve  Fit,  available  on 

the  Waterways  Experiment  Station's  (WES)  Honeywell  DPS8  computer  system. 

While  several  algebraic  equation  forms  fit  the  data  reasonably  well,  a  second 

order  polynomial  resulted  in  the  best  correlation  coefficient  for  both  rock 

types.  The  general  equation  form  may  be  expressed  as: 

2 


'uc 


a  +  b 


(-r^)  •  ■  {^) 


(15) 
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Corresponding  to  the  peak  deviator  stress  there  exists  a  unique  value  of 
axial  strain  (major  principal  strain)  such  that: 


£l  =  aL/L 


“^uc  ~  ^'"uc'^'^uc 


(13) 

(14) 


where 

ei  =  axial  strain  of  triaxial  compression  specimen  at  failure 

aL  =  total  axial  deformation  of  triaxial  compression  specimen  at 
failure 

L  =  original  length  of  triaxial  compression  specimen  at  failure 
=  axial  strain  of  unconfined  compression  specimen  at  failure 

aL  =  total  axial  deformation  of  unconfined  compression  specimen  at 
failure 

Luc  =  original  length  of  unconfined  compression  specimen  at  failure 

Note:  subscripts  denoting  failure  have  been  deleted  to  simplify 
presentation  of  terms. 

In  general  ei  >  For  the  unconfined  compression  case  where  03  =  0  the 

deviator  stress  a  .  =  oi  =  0 

Figure  5  shows  a  Mohr's  stress  circle  plot  intended  to  be  representative 
of  the  idealized  stress-strain  response  at  failure  shown  in  Figure  4.  The 
definition  of  the  pertinent  parameters  are  as  defined  in  Figure  5.  For  a 
given  rock  material  there  exists  a  unique  value  of  (unconfined  compression) 
and  ei  (triaxial  compression)  for  each  stress  circle.  In  developing  the 
proposed  empirical  strain  based  failure  criterion  presented  herein,  the  author 
assumed  that  the  failure  process  is  controlled  by  the  major  and  minor  princi¬ 
pal  stresses  and  o^,  and  that  the  intermediate  principal  02  has  no 
significant  influence  upon  this  process.  Although  this  assumption  over- 
simnlifies  actual  rock  failu'^e  be havior,  there  appears  to  be  sufficient  evidence, 
for  example  Jaeger  and  Cook  (1969),  to  suggest  that  the  influence  of  the  inter¬ 
mediate  principal  stress  can  be  ignored  without  introducing  unacceptably  large 
errors  for  practical  problems 


Triaxial  Test  Data 


The  author  analyzed  published  data  from  approximately  ninety  triaxial 
tests  on  intact  rock  specimens.  A  total  of  63  sandstone  and  23  granite  tests 
were  analyzed.  To  be  of  use  in  this  study  the  data  had  to  include  stress- 
strain  response  information.  This  information  usually  consisted  of  stress- 
strain  curves.  Data  were  taken  from  the  stress-strain  curves  reported  in 


Shear  Stress 


Tm  =  Maximum  shear  stress  =  (CTi-Qs)  /  2 

r  =  Shear  stress  along  failure  plane 

C3  =  Minor  principal  stress  (confining  stress) 

=  Major  principal  stress  (peak  axial  stress) 
o;,e  =  Unconfined  compressive  stress  (when  03  =  0, 

<Tl  -  0UC  ) 

On  =  Stress  normal  to  the  failure  plane 


Principal  Stress  Q 

Figure  5.  Idealized  Mohr's  stress  circles  and  failure 
envelope  for  tri axial  and  unconfined  compression  tests 
with  parameter  definitions 
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PART  III:  PROPOSED  STRAIN  BASED  FAILURE  CRITERION 


Background 

he  simplest  form  of  a  strain  based  failure  criterion  would  relate  the 
major  principal  strain  at  failure  to  the  principal  stresses  causing  failure 
for  a  given  intact  rock  type.  If  such  an  algebraic  expression  could  be  estab¬ 
lished  for  intact  rock,  then  with  suitable  case  history  comparisons  between 
predicted  and  observed  behavior,  the  expression  could  be  adjusted  for  rock 
mass.  This  part  of  the  report  will  examine  the  feasibility  of  developing  a 
strain  based  criterion  for  intact  rock. 

Specifically,  empirical  curves  best  defining  stress-strain  response  at 
peak  stress  will  be  fitted  to  tri axial  test  data  for  two  intact  rock  types  - 
sandstone  and  granite.  Intact  rock  was  chosen  for  this  study  because  of  the 
availability  of  tri axial  test  data  in  the  literature  and  because  failure 
mechanisms  for  intact  rock  are  simpler  than  mechanisms  associated  with  single 
plane  or  rock  mass  failure. 

Attributes 

At  the  onset  of  this  study  it  was  decided  that  any  proposed  empirical 
criterion  should  possess  four  basic  attributes: 

1.  The  criterion  should  be  simple  and  easy  to  understand. 

2.  Each  parameter  should  be  widely  accepted  by  the  rock  mechanics 
nractitioner, 

3.  Each  parameter  should  be  easy  to  measure  with  relatively  cheap 
and  repeatable  tests. 

4.  The  criterion  should  be  functional  in  that  all  parameters  are 
relevant  to  the  analyses  of  underground  excavation  behavior. 

'•'ith  these  attributes  in  mind  a  number  of  unsuccessful  correlation  attempts 
were  made  before  a  reasonable  correlation  was  obtained  between  stress  and 
strain.  The  scope  of  this  report  will  not  permit  a  discussion  of  the  unsuc¬ 
cessful  attempts.  Hence,  only  those  considerations  leading  to  a  reasonable 
correlation  will  be  discussed  herein. 

Parameter  Definitions 

Figures  4  and  5  define  the  necessary  parameters  used  in  the  development  of 
the  ''.train  based  criterion  proposed  herein.  Figure  4  shows  an  idealized 
stress-strain  response  for  simple  triaxial  compression  and  unconfined  com¬ 
pression  tests.  Failure  is  defined  by  peak  deviator  stress  (o^  =  (oi  -  03)  max). 
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Limitations.  The  empirical  rock  mass  failure  criterion  proposed  by  Hoek  and 
brown  (1980)  gives  an  approximate  method  for  estimating  the  strength  of  jointed 
rock  masses  in  terms  of  principal  stresses.  The  method  involves  estimating 
the  values  of  the  empirical  constants  "m"  and  "s"  from  a  description  of  the 
rock  mass  (summarized  in  Figure  3  and  Table  3).  These  estimates,  together 
with  an  estimate  of  the  uniaxial  compressive  strength  of  the  intact  pieces  of 
rock,  can  then  be  used  to  estimate  the  peak  deviator  stress  for  the 
jointed  rock  mass. 

Hoek  (1983)  states  that  -  "experience  in  using  the  values  (of  m  and  s) 
listed  in  Table  3  for  practical  engineering  design  suggests  that  they  are 
somewhat  conservative."  In  addition  to  being  somewhat  conservative,  the 
criterion  is  valid  only  for  brittle  failure.  Schwartz  (1964)  indicated  that 
for  Indiana  Limestone  a  transition  from  brittle  to  ductile  behaviour  appeared 
to  occur  at  a  principal  stress  ratio  of  approximately  01/03  =  4.3.  Mogi 
(1966)  concluded  that  the  transition  for  most  intact  rock  occurs  at  an  average 
principal  stress  ratio  01/03  =  3.'".  Considering  the  curve  fitting  procedures 
and  the  variety  of  interpretations  employed  in  defining  the  transition  Hoek 
(1983)  suggested  a  rough  conservative  rule-of-thumb  to  define  the  upper  limits 
of  stress  for  brittle  failure  of  intact  rock.  That  rule  states  that  the  minor 
principal  (confining)  stress  must  always  be  less  than  the  uniaxial  compressive 
strength  of  the  material.  While  suitable  guidance  for  rock  mass  behavior  is 
not  available  it  seems  reasonable  to  conclude  that  brittle-ductile  transitions 
for  a  rock  mass  would  be  of  the  same  order  of  magnitude  as  for  intact  rock. 
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GEOMECHANICS  CLASSIFICATION  (RMR) 

3.  Criterion  for  in  situ  strength  of  rock 
(after  Hoek  and  Brown  1980) 


form  of  the  criterion  applicable  to  the  proposed  constitutive  model.  Hence, 
single  discontinuities  will  not  be  addressed  herein. 

"m"  and  "s"  for  rock  mass.  The  presence  of  one  or  more  sets  of  discontinuities 
in  a  rock  mass  will  cause  a  reduction  in  the  values  of  both  "m"  and  "s". 

Unlike  intact  rock,  relatively  few  sets  of  reliable  triaxial  test  data  for 
jointed  rock  are  available.  As  such,  the  choice  of  "m"  and  "s"  for  a  given 
rock  mass  must  be  based  on  the  limited  data  available  as  well  as  back-analysis 
of  documented  cases  of  rock  mass  failures  (as  they  become  available)  and  a  good 
measure  of  judgement.  The  only  known  set  of  triaxial  test  data  available  for 
undisturbed  naturally  jointed  rock  was  that  obtained  by  Jaeger  (1970)  for  the 
Panai'r.?  andesite  from  the  Bougainville  open  pit  mine,  Papua-New  Guinea. 

Figure  3  shows  plots  on  the  constant  "s"  and  the  ratio  of  "m/m^,"  in 
which  "m^"  is  the  value  of  "m"  for  intact  rock  material,  against  the  Rock  Mass 
Rating  (RMR)  system  developed  by  Bieniawski  (1973)  and  the  Q-System  developed 
by  Barton,  Lien  and  Lunde  (1974).  The  RMR  and  Q  classification  ratings  are 
estimates  for  the  various  categories  of  the  Panguna  andesite  listed  in  Table 
2.  An  inspection  of  Table  2  indicates  that  5  data  points  for  each  classifi¬ 
cation  rating  are  available  for  "m"  with  none  above  an  RMR  of  approximately  46. 
One  point  for  each  classification  rating  (at  an  RMR  of  approximately  46)  is 
available  for  "s"  (plus  s  =  1.0  for  intact  rock).  Straight  lines  were  drawn  on 
Figure  3  to  give  approximate  relationships  between  "s"  and  "m/m^"  and  the 
classification  rating.  Table  3  lists  approximate  relationships  between  rock 
mass  quality  and  the  material  constants  suggest  by  Hoek  and  Brown  (1980)  and 
Hoek  (1983). 

Priest  and  Brown  (1983)  recently  proposed  a  probabilistic  method  for 
analysing  the  stability  of  slopes  excavated  in  fractured,  variable  rock 
masses.  The  method  is  based  on  a  Monte  Carlo  simulation  routine  with  the  use 
of  the  Janbu  method  of  slices  for  which  strength  parameters  are  derived  by  the 
Hoek  and  Brown  failure  criterion.  Priest  and  Brown  (1983)  used  mathematical 
expressions  to  determine  values  of  the  material  constants  "m"  and  "s"  from 
Figure  3  as  follows; 


R  =  exp 


Rf^R  -  95 
13.4 


m  =  R  m-j 
s  =  exp 


RMR  -  100 
6.3 


Uhere  all  parameters  are  defined  above. 


(10) 

(11) 

(12) 


nonlinear  failure  envelope  predicted  by  classical  Griffith  crack  theory 
for  plane  compression  (Hoek  1968)  and  by  using  a  process  of  trial  and  error. 
The  failure  criterion  is  expressed  by  the  following  empirical  relationship 
between  the  principal  stresses  at  failure: 

^%c  ""  ”3  /^uc 

Where 

oi  =  the  major  principal  stress  at  failure 
03  =  the  minor  principal  stress  at  failure 

=  the  uniaxial  compressive  strength  of  the  intact  rock  material 

m  &  s  =  material  constants  that  depend  on  the  properties  of  the 
rock  and  on  the  extent  to  which  it  had  been  broken  before 
being  subjected  to  the  failure  stress  oi  and  03. 

Equation  7  can  be  rewritten  as  follows: 

oi  =  03  +  (m  ajjj.o3  +  (8) 

The  scope  of  this  report  will  not  permit  a  detailed  discussion  of  the 
criterion  development.  For  details,  the  reader  is  referred  to  Hoek  and  Brown 
(1980)  and  Hoek  (1983).  However,  for  continuity  of  following  discussions 
material  constants  "m"  and  "s"  will  be  briefly  discussed. 

"m”  and  "s"  for  intact  rock.  For  isotropic  intact  rock,  the  material  constant 
"s"  is  defined  as  equal  to  1.0.  Hence,  for  intact  rock  Equation  8  becomes: 

oi  =03+  (ma^j^  03  +1.0 

A  regression  analysis  was  then  used  to  determine  the  value  "m"  giving  the  best 
fit  of  Equation  9  to  sets  of  normalized  triaxial  data  (e.g.  plots  of 
oi  /  versus  aj  /  o^^)  for  a  particular  rock  type.  Hoek  (1983)  summarizes 
the  best  fit  procedure.  The  value  of  "m"  for  a  given  rock  type  is  dependent 
on  the  crystal  fabric  of  the  rock  material  and  varies  from  approximately  5  to 
30  (Hoek  and  Brown  1980). 

"m”  and  "s"  for  discontinuous  rock.  The  material  constants  "m"  and  "s"  for 
single  discontinuous  planes  are  not  constant  as  they  are  for  intact  rock,  but 
vary  with  the  orientation  of  the  plane  of  weakness.  It  is  unlikely  that  the 
criterion  would  be  particularly  useful  in  practice  because  of  the  somewhat 
complicated  procedures  required  to  define  "m"  and  "s"  nor  is  the  single  plane 


cannot  be  expected  to  be  representative  of  modulus  defined  at  peak  stress 
(i.e.  slope  of  line  OB  in  Figure  2).  Finally,  the  assignment  of  a  unique 
value  for  the  deformation  modulus  ignores  the  fact  that  modulus  varies  with 
confining  stress  (minor  principal  stress)  as  indicated  in  Figure  1. 

Proposed  Constitutive  Model 


Minimum  Requirements 

A  completely  general  and  ideal  constitutive  model  would  consist  of  an 
algebraic  expression,  or  set  of  algebraic  expressions,  that  defines  the  complete 
rock  mass  stress-strain  curve  (to  include  prefailure,  failure  and  postfailure). 

As  a  minimum,  a  model  convenient  for  underground  excavations  would  need  to  be  a 
function  of  the  major  principal  strain,  the  major  and  minor  principal  stress,  and 
the  rock  mass  conditions.  Development  of  such  a  model  would  be  a  formidable, 
if  not  impossible,  task.  Nevertheless,  the  practicing  geotechnical  engineer 
is  constantly  searching  for  improved  techniques  for  estimating  stress-strain 
response. 

Previous  discussions  on  failure  criteria  indicated  that  a  number  of 


criteria  are  available  for  estimating  stress  at  failure.  However,  no  suitable 
criteria  are  available  for  estimating  strain  at  failure.  The  development  of  a 
strain  based  failure  criterion  would  provide  a  set  of  algebraic  expressions 
for  estimating  deformation  modulus  up  to  failure  as  defined  by  peak  deviator 
stress  and  strain  at  failure  (i.e.  Equation  4).  The  general  form  of  expres¬ 
sions  for  peak  deviator  stress  and  strain  at  failure  would  be  as  follows; 


f(oi,  03,  rock  mass  conditions) 
f(ai,  03,  rock  mass  conditions) 


(5) 

(6) 


While  such  a  formulation  (Equation  4)  would  not  be  ideal  in  that,  if  required. 


postfailure  stress-strain  behavior  would  still  have  to  be  defined  by  some 
other  model,  it  would  provide  a  significant  advancement  for  estimating  the 
initial  moduli  values. 


'Jhile  no  suitable  strain  based  criteria  existed  for  rock  materials,  a 


stress  based  criterion  applicable  to  intact  rock,  discontinuous  rock  and  rock 


mass  is  available.  This  stress  based  criterion,  proposed  by  Hoek  asiu  Brcwn 
(1980),  meets  the  minimum  requirements  as  specified  in  Equation  5. 


Hoek  and  Brown  Failure  Criterion 

The  Hoek  and  Brown  failure  criterion  was  developed  by  analogy  with  the 
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€/€«:=  0.26  0.73 +  0.02  (5^) 
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A 
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■ 

Westerly  Granite 
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Figure  7.  Plot  of  normalized  axial  strain  versus  normalized 
maximum  shear  stress  for  granite 


Correlations  Between  Minor  Principal 
Stress  (Confining  Stress)  and  Axial  Strain 


Curve  Fittinc 


Curve  fitting  procedures  for  the  normalized  minor  principal  stress 
(03/0  )  and  normalized  axial  strain  (e/e  )  data  were  similar  to  those 

described  above.  The  data  are  given  in  Table  A1  of  Appendix  A.  The  Curve 
■^it  program  analysis  resulted  in  a  best  fit  second  order  polynomial  of  the 
general  form: 


+  b  —  +  c  — \ 

%c  \  I'C  j 


For  the  unconfined  compression  test,  where  03=0;  ,  Equation  17  is  valid. 
Hence,  from  Equations  17  and  21,  the  minimum  value  of  the  normalized  axial 
strain  parameter  is  as  follows: 


_  uc  _ 


a  =  1.0 


Values  of  the  constants  b  and  c  in  Equation  21  were  then  adjusted  by  trial  and 
error  to  provide  a  visual  best  fit  of  the  data  and  satisfy  the  minimum  value 
criterion  defined  in  Equation  22. 

Sandstone  Data 

'^igure  8  shows  a  plot  of  normalized  axial  strain,  ,  vs  normalized 

minor  principal  stress,  03/0,^^  ,  for  the  sandstone  data  given  in  Table  A1  of 
Appendix  A.  The  expression  for  the  line  best  ftting  all  the  data  is  as 
follows: 


=  1.00  + 


(1.20.  +  0.25  03/0^^^) 


Granite  Data 

'"igure  9  shows  a  plot  of  normalized  axial  strain,  ,  vs  normalized 

minor  principal  stress,  03/0^^  ,  for  the  granite  data  given  in  Table  A1  of 
Appendix  A.  The  expression  for  the  line  best  fitting  all  the  data  is  as 
follows: 


=  1.00 


+  -5^  /4.3  -  0.3 
%c  ^  °uc  / 
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Confining  Stress/Unconfined  Compressive  Strength  —  03/oue 

Figure  8.  Plot  of  normalized  axial  strain  versus  normalized 
minor  principal  stress  for  sandstone 
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Deformation  Modulus 


The  deformation  modulus  as  defined  by  Equation  4  can  be  rewritten  as 
follows: 

Ed=  (25) 

Hence,  in  order  to  predict  deformation  modulus  one  must  be  able  to  predict  the 
principal  stresses  at  failure  and  axial  strain  at  failure.  Hoek  and  Brown's 
(1980)  failure  criterion,  as  defined  in  Equation  7,  offers  a  viable  means  of 
predicting  the  normalized  principal  stress  difference  at  failure.  For  con¬ 
venience  Equation  7  can  be  expressed  as: 

=  (m  03  t  s)’/2  ^26) 

Strain  based  criterion  proposed  herein  allows  estimates  of  strain  at 
failure  for  intact  sandstone  and  granite.  As  proposed,  normalized  axial  strain 
at  failure  is  expressed  as  functions  of  either  normalized  maximum  shear  stress 
or  normalized  minor  principal  stress.  Strain  expressions  in  terms  of  the 
minor  principal  stress  renresent  the  more  convenient  form  for  two  reasons. 
First,  note  that  the  Hoek  and  Brown  criterion.  Equation  26,  expresses  the 
normalized  deviator  stress  as  a  function  of  the  normalized  minor  principal 
stress  (as  well  as  the  constants  m  and  s). 

Expressing  the  normalized  strain  as  a  function  of  normalized  minor  princi¬ 
pal  stress  would  reduce  the  number  of  stress  parameters  in  Equation  25  to  one; 
that  being  the  minor  principal  stress.  Finally,  in  underground  excavations, 
the  minor  principal  stress  is  known  or  assumed  to  within  reasonable  tolerances 
at  two  boundaries;  at  the  excavation  surface  and  at  some  unknown  distance 
within  the  wall  rock  unaffected  by  the  relaxation  process.  The  major  prin¬ 
cipal  stress,  however,  is  seldom  known  except  at  some  unspecified  distance 
within  the  wall  rock.  Numerical  analyses  require  a  known  boundary  condition, 
such  as  at  the  excavation  surface  where  both  the  location  of  the  boundary  and 
minor  principal  stresses  are  known.  Hence,  deformation  moduli  expressed  in 
terms  of  the  minor  principal  stress  are  more  conducive  to  numerical  modeling. 

A  dimensionless  form  of  the  deformation  moduli  expressed  as  a  function  of 
the  minor  principal  stress  can  be  obtained  for  intact  sandstone  or  granite  by 

substituting  the  left  hand  sides  of  Equation  26  and  Equations  23  or  24  (for 
sandstone  or  granite  respectively)  into  Equation  25  as  follows: 
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(27) 


rewritten  in  terms  of  the  normalized  minor  principal  stress  ratio  the 
equation  for  sandstone  becomes: 


where  m  =  15  (from  Table  3) 
and  for  granite 

Sc 

°uc  1.0  +  ^  (4.3-0. 3  o^) 


(29) 


where  m  =  25  (from  Table  3) 

Equations  28  and  29  are  plotted  in  Figure  10.  It  is  interesting  to  note  that 
for  small  increases  in  the  03/0^^^  ratio  both  rock  types  exhibit  sharp  in¬ 
creases  in  the  dimensionless  ^uc'^'^uc  ®  clearly  defined  peak. 

Granite  peaks  at  a  smaller  03/0^^.  ratio  (03/0^^.  =  0.2)  than  sandstone 
(03/Oyc  ~  0.5).  With  increasing  03/°^.  ratio  past  peak  the  E^j  ratio 

decreases.  However,  note  that  for  03/0^^  ratios  greater  than  approximately 
1.0  Equation  29  for  intact  granite  predicts  a  lower  modulus  that  would  be 
obtained  from  unconfined  compression  tests  (i.e.  for  unconfined  compression 
Ed  =  1.0).  Sandstone,  on  the  other  hand,  requires  a  03/0^^^,  of  greater 

than  approximately  3.0  before  the  E^  '^uc^^uc  falls  below  1.0. 
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Figure  10.  Plot  of  normalized  deformation  modulus  versus 
normalized  confining  stress  for  sandstone  and  granite 
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PART  IV:  SUMMARY  AND  CONCLUSIONS 
Summary 

The  constitutive  model  proposed  in  this  reoort  equates  the  deformation 
modulus  to  the  ratio  of  the  deviator  stress  vs  axial  strain  at  failure.  For 
the  model  to  be  of  practical  significance  requires  failure  criteria  which 
define  both  the  deviator  stress  at  failure  and  strain  at  failure.  The  Hoek 
and  Brown  (1980)  failure  criterion  offers  a  potentially  powerful  tool  for 
defining  the  deviator  stress  in  that  the  criterion  expresses  the  deviator 
stress  as  a  function  of  minor  principal  stress  normalized  with  respect  to  the 
unconfined  compressive  strength  and  material  constants  "m"  and  "s."  The 
material  constants  vary  depending  upon  rock  conditions  ranging  from  intact  to 
a  highly  jointed  rock  mass.  Hoek  and  Brown  (1980)  offer  the  only  criterion 
that  addresses  rock  mass  failure. 

A  detailed  literature  search  failed  to  locate  a  criterion  suitable  for 
rock  materials  that  defines  strain  at  failure.  In  an  effort  to  fill  this  void 
a  feasibility  study  was  conducted  to  evaluate  the  possibility  of  developing  a 
strain  based  criterion  for  rock  materials. 

Stress-strain  responses  for  two  rock  types,  sandstone  and  granite,  were 
evaluated.  Reasonable  correlations  were  found  to  exist  between  axial  strain 
at  failure  normalized  with  respect  to  unconfined  compressive  axial  strain  at 
failure  and  the  following: 

1.  Maximum  shear  stress  normalized  with  respect  to  the  unconfined 
compressive  strength. 

2.  Minor  principal  stress  (confining  stress)  normalized  with 
respect  to  the  unconfined  compressive  strength. 

The  Hoek  and  Brown  criterion  and  the  latter  correlation  provides  an  expression 

for  estimating  deformation  modulus  as  a  function  of  the  minor  principal  stress. 

Conclusions 


This  study  has  demonstrated  the  feasibility  of  developing  empirical 
relationships  defining  strain  at  failure  for  rock  materials,  at  least  for 
intact  sandstone  and  granite.  Practical  application  requires  that  the  rela¬ 
tionships  be  extended  to  rock  mass  conditions.  Extension  of  the  criteria 
can  be  accomplished  in  two  ways.  The  obvious  approach  would  consist  of  a 
series  of  triaxial  tests  on  specimens  of  sufficient  size  to  be  representative 
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of  various  rock  mass  conditions.  Such  an  approach  is  prohibitively  expensive 
and  would  require  large  diameter  high  load  triaxial  equipment  not  commonly 
available.  A  second  alternative  approach  would  consist  of  using  numerical 
analysis  of  existing  case  histories  to  adjust  the  algebraic  form  of  criteria 
so  that  predicted  behavior  matches  observed  behavior.  Such  an  approach 
would  unavoidably  impose  an  additional  assumption,  that  assumption  being  the 
validity  of  the  numerical  technique. 

Work  is  in  progress  at  the  Waterways  Experiment  Station  in  an  effort 
to  develop  a  strain  base  failure  criterion  applicable  to  rock  mass  conditions. 
If  successful,  the  criterion  will  significantly  increase  our  ability  to  pro¬ 
vide  meaningful  constitutive  relationships  between  stress  and  strain  necessary 
for  rational  design  of  underground  excavations. 
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Table  1 

Major  Failure  Criteria  for  Rock  Materials 
(Bieniawski  1984) 


Name  of  Criterion 

Proposer 

Date 

Reference 

Coulomb-Navi er 

Coulomb 

1773 

Jaeger  and  Cook  1979 

Maximum  shear  stress 

Tresca 

1864 

Nadai  1950 

Maximum  principal  stress 

Rankine 

1869 

Nadai  1950 

Maximum  elastic  strain 

St.  Venant 

1870 

Nadai  1950 

Constant  elastic  energy 
of  deformation 

Bel trami 

1885 

Nadai  1950 

Shear  failure 

Mohr 

1900 

Jaeger  and  Cook  1979 

Constant  octahedral  shearing 
stress 

Huber 

Hencky 

1904 

1920 

Jaeger  and  Cook  1979 

Second  invariant  of  stress 
deviation 

von  Mises 

1913 

Jaeger  and  Cook  1979 

Griffith-original 

Griffith 

1921 

Griffith  1921  ,1924 

Statistical  failure  theory 

Wei  bull 

1939 

Weibull  1952 

Fracture  toughness 

Irwin 

1960 

Irwin  1960 

Griffith -modified 

McClintock-Walsh 

1962 

McClintock  and  Walsh 

1962 

Hoek  and  Bieniawski 

1965  I 

Griffith-extended 

Murrel 1 

1963 

Murrell  1963 

Empirical  rock  mass  strength 

Hoek-Brown 

1980 

Hoek  and  Brown  1980 

NGi  classification  csir  classification 


The  Rock  Mass  Rating  and  Q-System  Classifications  of 
Panguna  Andesite  Samples  (Hoek  and  Brown  1980) 
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0 

0.0006 

0 

Intact  strength  rating 

15 

15 
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15 

15 

15 

CI  RQO  rat  i ng 

20 
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3 

3 

3 

3 

Joint  spacing  rating 
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5 

u.  Joint  condition  rating 
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2  Groundwater  rating 
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to 
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Joint  let  number  Jp 
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9 

t2 
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18 

20 

Joint  roughness  nuniber  J^ 

1* 

3 

1.5 

1 

1 

1 

Zt  Joint  alteration  number  J^ 

0.75 

2 

1) 

6 

8 

2  Joint  water  reduction  Jy, 

1 

I 

t 

t 

1 

1 

—  Stress  reduction  factor  5RF 

I 

2.5 

5 

7.5 

to 

10 

^  NGl  quality  index  Q 

533 

0.67 

0.06 

0.02 

0.009 

0.006 

Table  3 

Approximate  Relationship  Between  Rock  Mass  Quality, 
Rock  Type,  and  Material  Constants  (Hoek  1983) 
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tXTACT  ROCT  SAMPLES 

laboratory  alee  aaxplca  free 
froa  pre*eaiaclat  fracturca 

■icniavaki  (CSIR)  rating  100 
Barton  (NCI)  racing  SOO 

VERY  GOOD  QUALITY  ROCR  HASS 

Tightly  interlocking  undiatur- 
bed  rock  with  rough  unveachered 
jolnca  apaccd  at  I  to  3  a 

Bienlavakl  <CS1R)  racing  65 
Barton  (HCI)  rating  lOQ 

GOOD  quality  ROCR  MASS 

Preah  to  allghcly  weathered 
rock,  slightly  diaturbcd  with 
jointa  apaced  at  1  to  3  a 

Bienlawaki  (CSIR)  rating  6S 
Barton  (MCI)  rating  10 


0.004  a  -  0.004 


PAIR  QUALITY  ROCX  MASS 

Several  aeca  of  aodcratcly 
vcatbered  Jointa  apaced  at 
0. 3  to  1  ■  •  diaturbed 

Bienlawaki  (CSIR)  rating  44 
Barton  (NCI)  rating  I 

POOR  QUALITY  BOCK  MASS 

Huaeroua  leeatbered  jointa  at 
30  to  SOOoa  with  aoac  gouge. 
Clean,  coapacted  rockfili 

Bienlawaki  (CSIR)  rating  23 
Barton  (NCI)  rating  0.1 


a  •  0  0001  a  "  0.0001  a  •  0.0001  a 


a  •  0.34  a  •  0.50 
I  -  O.DOOl  a  •  0.0001 


0.  OOOOl  1  a  •  O.OOC'Cl  a  -  0. 0000 1  a  •  0. 0000 1  la  -  0.00001 


VERT  POOR  QUALITY  ROCK  MASS 

Ihiacroua  heavily  weathered 
Jointa  apaced  at  50aai  with 
gouge.  Vaatc  rock 


Bienlawaki  (CSIR)  racing  3 
Barton  (HCX)  racing  0.01 
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